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EXECUTIVE  SIMIARY 


The  Wideband  Satellite  data  presently  constitute  the  most  extensive 
collection  of  global  Intensity  and  phase  scintillation  data  available. 

The  data  are  being  used  to  evaluate  the  scintillation* Induced  performance 
degradation  In  a variety  of  systems.  For  a low  orbiting  satellite  such 
as  Wldebax.a,  however,  the  structure  of  the  scintillation  data  Is 
critically  dependent  on  the  propagation  geometry.  Thus,  to  use  the 
Wideband  satellite  data  effectively,  we  must  have  a model  to  separate 
the  purely  geometrical  effect^'  from  the  true  source  variations. 

In  this  report  we  show  that  an.  appropriately  located  equivalent 
phase-changing  screen  accurately  reproduces  the  measured  phase  and 
Intensity  scintillation  from  Individual  Wideband  passes.  We  assume  that 
the  three-dimensional  spectral  density  of  the  Ionospheric  Irregularities 
has  the  form  '.  The  spectral  Index  parameter  v and  the  strength 

of  turbulence  are  determined  by  using  the  model  calculations  and  the 
suDBoary  parameters  that  are  routinely  obtained  from  the  phase  scintillation 
data  by  fitting  a nurve  of  the  form  Tf  ^ to  the  measured  phase  spectral 
density  fiinctlon. 

The  measured  average  p Index  values  fall  In  the  range  2 < p < 3; 
however,  the  data  show  a systematic  difference  between  the  auroral  zone 
and  the  equator.  For  the  auroral-zone  data,  the  p Indices  are  generally 
smaller  and  more  variable  than  the  p Indices  for  the  equatorial  data, 
which  give  values  closer  to  the  nominal  value  p =•  3.  Nonetheless,  the 
KWajaleln  dots  give  somewhat  smaller  p values  them  the  Ancon  data. 

The  free  parameters  In  the  phase- screen  model  are  the  height  of  the 
phase  screen  emd  the  Irregularity  axial  ratios  along  emd  transverse  to 
the  magnetic  field.  We  have  assumed  that  the  satellite  scan  velocity, 
which  Is  known,  dominates  the  irregularity  drift  component.  Self- 
consistent  fits  to  the  measured  rms  phase  and  generally  require  the 
equivalent  phase  screen  to  be  in  the  F region  for  both  the  auroral  and 
equatorial  data. 


The  equatorial  data  are  conaiatant  with  rod- like  Irregularities  with 
a mlnlmuat  axial  ratio  of  20:1.  The  auroral-zone  data  are  generally 
consistent  with  sheet-llke  structures  aligned  along  L-shells  at  least  for 
propagation  paths  wi.thln  the  auroral  oval.  A detailed  analysis  of  the 
Irregularity  anisotropy  Is  being  Independently  pursued  by  using  the 
Wideband  spaced  recislver  data. 

The  phasA- screen  formulas  have  been  greatly  simplified  by  taking 
limits  as  the  outer  scale  approaches  Infinity  and  the  Inner  scale 
approaches  zero.  The  Justification  for  using  these  approximations  Is 
that  no  evidence  of  systematic  departures  from  the  power- law  spectral 
form  has  been  found  in  either  the  Wideband  data  or  any  other  reported 
phase  scintillation  data.  Thus,  the  cutoff  scales  are  outside  the  range 

of  the  scale  sizes  t:hat  affect  the  data,  and  they  are  properly  excluded 
from  the  analysis. 

We  have  also  used  the  estimated  turbulence  levels  to  determine  the 
rms  electron  density  perturbation  that  would  be  measured  by  an  In- situ 
probe.  We  have  found,  for  example,  that  significant  gigahertz  scintil- 
lation can  be  accounf.ed  for  by  purturbatlons  with  rms  electron  density 
levels  between  10  cl/m  and  10^^  eiVm^  over  a 200-km  propagation  path. 

Smaller  levels  distributed  over  a larjuer  path  would,  of  course,  produce 
the  same  result. 
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I XWTROOUCTION 


Phapd  changing  acraenb  have  been  uaau  to  model  optical  and  iadlo* 
wave  scintillation  phenomena  throughout  the  development  of  the  theory  of 
scintillation  (Booker  at  al.^  1950;  Bramley,  19S4;  Ratcllffe,  1956; 
Salpeter,  1967;  Cronyn,  1970;  Rufenach,  1975;  Rumsey,  1975;  B».ickley, 
1975).  The  phase-screen  model  Is  appealing  because  of  Ita  conceptual 
and  nnalytlcal  simplicity.  It  Is  often  argued,  hovrever,  that  the  phase- 
screen  model  cannot  Accnrajeiy  describe  scintillation  data  because  non- 
negllglble  amplitude  variations  Invariably  develop  within  the  scattering 
medium  (see  for  example  Wernlk,  1976;  Crane,  1977). 

While  such  arguments  are  eseentuslly  correct,  one  feels  Intuitively 
that  an  appropriately  located  equivalent  phase  screen  should  accurately 
reproduce  the  average  properties  of  scintillations  that  develop  both 
within  and  beyond  an  extended  medium.  Indeed,  Bramley  (1977)  has 
recently  demonstrated  by  direct  computation  that  such  an  equivalent 
phase  screen  Is  surprisingly  accurate.  For  Irregularities  with  a power- 
law  spectral  distribution,  Bramley' s calculations  show  that  the  errors 
made  in  computing  the  amplitude  ana  phase  scintillation  levels  by  using 
an  equivalent  phase  screen  amount  to  only  a f( percent  even  If  the 
receiver  lies  within  the  scattering  medium. 

In  this  report  we  shall  demonstrate  that  the  phase-sezeen  model 
accurately  reproduces  both  the  level  and  the  propagation  geometry  depen- 
dence of  Ionospheric  scintillation  data  obtained  from  the  Wideband  satel- 
lite. Data  have  been  acquired  from  auroral  zone  and  equatorial  stations. 
The  auroral-zone  data  present  a complicated,  rapidly  changing  propagation 
geometry  combined  with  extreme  source  variations.  Ihe  equatorial  data 
provided  excmplcs  of  the  gigahertz  scln:  ..xlatlon  phenomena  that  have 
recently  received  considerable  attention  (Basv  et  al.,  1976;  Costa  and 
Kelley,  1976). 


We  ehall  ccn.vlder  only  weak  scattering  In  this  report.  This  is  not 
overly  restrictive  because  saturated  scintillation  at  L-band  Is  observed 
only  occasionally.  Nonetheless,  even  under  conditions  of  strong  scattering, 
the  phase- screen  model  remains  the  simplest  means  of  obtaining  flexible 
yet  tractable  results,  particularly  in  a power- law  scattering  environment 
(Gochelashvlli  and  Shlshov,  1971;  Taylor,  1972;  and  Rumsey,  1975). 

Phase  scintillation  data  are  used  in  the  model  calculations 
essentially  to  determine  the  rms  phase  le'  1 for  the  phase- changing 
screen.  Indeed,  the  potential  of  differential  phase  measurements  has 
been  recognized  for  some  time  as  an  accurate  means  of  measuring  relative 
and  absolute  tot  il  (integrated)  electron  content  (TEC)  (Leitlnger  et  al., 
1975). 

The  principle  of  the  differential  phase  technique  is  simply  that 
large-scale  phase  variations  are  directly  proportional  to  wavelength 
times  the  TEC  along  the  propagation  path.  The  more  rapidly  changing 
phase  scintillations  can  be  thought  of  as  TEC  noise.  To  some  extent, 
diffraction  effects  r.  jst  ultimately  upset  the  simple  proportionality 
between  phase  and  wavelength  times  TEC.  An  Important  finding  from  the 
Wideband  data  is  that  under  conditions  of  weak  scattering,  diffraction 
effects  in  the  phase  data  are  .negligible  to  time  scales  substantially 
shorter  than  the  characteristic  intensity  fading  period. 

If  diffraction  effects  were  indeed  Important,  one  would  expect  the 
phase  scintillation  to  be  sensitive  to  changes  in  Fresnel  radius  and  thus 
vary  r jnlinearly  with  wavelength.  Yet,  we  have  consistently  observed  a 
nearly  linear  wavelength  dependence  of  rms  phase  except  under  the  most 
severe  scattering  conditions  (Fremouw  et  al.,  1978).  We  conclude  that 
under  conditions  of  weak  scattering  the  diffraction  effects  in  the  phase 
data  can  be  largely  ignored.  We  shall  demonstrate  this  in  Section  II  by 
comparing  simultaneously  recorded  phase  data  at  different  frequencies. 

A second  important  finding  in  both  the  Wideband  data  and  other 
differential  phase  measurements  (Crane,  1976)  is  that  the  differential 
phase  spectrum  admits  a power- law  representation  with  no  intrinsic  low- 
frequency  cutoff  (outer  scale).  The  range  of  the  power-law  continuum  has 
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recently  been  discussed  by  Yeh  and  Liu  (1977).  The  fact  that  we  are 
deall'ig  with  a power- law  contlnutim  with  an  Inaccessible  outer  scale  has 
Important  consequences  for  data  Interpretation.  Indeed,  neither  an 
Inner  scale  nor  an  outer-scale  cutoff  has  been  detected  In  any  scintil- 
lation data  reported  to  date.  It  follows  that  such  parameters  need  not 
enter  a properly  formulated  scintillation  theory. 

Tlie  Inner- scale  cutoff  does  not  affect  the  amplitude  or  phase 
scintillation  data  because  It  occurs  well  below  the  nolje  level.  In 
the  Appendix  we  present  a simple  formula  that  allows  us  to  correct  the 
Intensity  scintillation  index  for  noise  contamination.  Nc  correction 
is  necessary  for  the  phase  data. 

If  the  outer-scale  cutoff  Is  large  compared  to  the  Fresnel  radius, 
the  diffraction  process  removes  any  influence  that  large  scale  structure 
might  have  in  the  Intensity  scintillation  data.  In  the  phase  data  the 
detrending  procedure  that  is  used  to  remove  TEC- induced  slow  phase  trends 
dictates  the  low-frequency  cutoff  of  the  spectrum.  Under  conditions  of 
both  weak  and  strong  scattering,  Fresnel  filtering  introduces  an  intrinsic 
low-frequency  cutoff  in  intensity  sc  inti  Up*:  ion  data.  Thus,  since  the 
outer  scale  has  not  been  detected  in  phase  scintillation  data,  it  follows 
a fortiori  that  the  outer  scale  has  no  direct  Influence  on  intensity 
scintillation  data.  We  shall  see  in  Section  III  that  if  one  ttdees  this 
fact  into  account  at  the  outset,  the  theory  is  considerably  simplified. 

In  Section  IV  we  apply  the  power- law  phase- screen  theory  to  Wideband 
equatorial  and  auroral  data.  We  show  that  whenever  the  model  accurately 
reproduces  the  rms  phase  data,  there  is  a set  of  anisotropy  parameters 
and  an  equivalent  phase  screen  height  that  will  accurately  reproduce  the 
measured  S^  values  as  long  as  0.4  For  larger  values  the  weak- 

scatter  theory  overestimates  S^. 

By  using  the  phase- screen  model  we  can  estimate  the  in- situ  rms 
electron  density  perturbation  level  that  is  consistent  with  the  scintil- 
lation data.  In  Section  V '.‘w  show  that  gigahertz  scintillation  can  be 
easily  explained  with  rms  electron  density  perturbations  between 
10  ^ el/m  and  10  el/m  distributed  over  a 200- km  layer. 
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II  PEASE  SCINTILLATION 

In  this  section  we  shall  first  review  a general  model  that  completely 
characterizes  the  structure  of  ths  phase  scintillation  exclusive  of 
diffraction  effects.  We  shall  then  compute  the  form  of  the  one- dimensional 
phase  power  spectrum,  which  can  actually  be  measured.  Prom  the  form  of 
the  phase  power  spectrum  we  deduce  the  relationship  that  converts  the 
measured  power  In  a given  frequency  Interval  to  an  electron  der^'ity 
perturbation  level  In  the  corresponding  spatial  frequency  regime. 

We  also  calculate  the  fom  of  the  sieasurable  rms  phase  and  discuss 
the  consequences  of  the  Inaccessible  outer  scale  wavenumber.  Finally,  we 
present  some  examples  that  show  the  extent  to  which  diffraction  effects 
are  detectable  In  phase  data  under  conditions  of  weak  scattering. 

The  zeroth-order  approximation  to  the  differential  phase  60  is  given 
by  the  Integral  along  the  propagation  path 


■r  X 1 ^ dt 
el  e 


(l  - (f/f  )^)  + terms  that  depend 
^ ^ on  Xz 


where  Is  the  classical  electron  radius,  X Is  the  wavelength  (fX  » c, 
where  c Is  the  velocity  of  light),  f^.  Is  the  reference  frequency,  and 
/^sNg  Is  the  local  electron  density  perturbation.  For  the  moment,  let  us 
assume  that  f^  is  Infinite  and  that  Xz-dependent  terms  In  Eq.  (1)  are 
Indeed  negligible. 

In  Rlno  and  Fremouw  (1977)  It  Is  shown  that  the  phase  autocorrelation 
function  derived  from  Eq.  (1)  has  the  general  form 


^60  ^^s^  “ 


■tan  0 


. K)  cos  (IC  • ^ ) 


where  L is  the  layer  thickness 


tan  0 


•nd 

■ (co«  ^p,  sin  (p)  . (4) 

The  angles  9 and  9 are,  respectively,  the  senlth  and  aiagnetlc  aclnuth 
angles  of  the  propagation  vector  k.  Ws  note  that  lies  along  the 

-4  ** 

horlsontal  projection  of  k.  The  z axis  of  the  reference  coordinate 
system  Is  downward-dlr<'icted  and  the  xz  plane  contains  the  local  geo- 
magnetic meridian. 

We  note  that  Eq.  (2)  Is  a fully  three-dimensional  characterization 

of  the  phase  structure.  The  £iZ  dependence  Is  contained  in  In 

-•  * 

Eq.  (2)  *^^(*^**^2^  three-dimensional  spectral-density  function 

(SDF)  of  the  Irregularities.  In  Rino  and  Fremouw  (1977),  It  Is  shown 

that  for  a fairly  general  anisotropy  model,  the  SDF  In  Eq.  (2)  has  the 

form 


-tan  9 


K) 


ah  {Mb  Q (AkJ  + BK  It  + ckJ)  . 
e X X y y 


(5) 


The  parameters  a and  b are  axial  ratios  along  and  transverse  to  the 
principal  Irregularity  axis.  The  coefficients  A,  B,  and  C depend  on  the 
propagation  angles  relative  to  the  principal  irregularity  axis  (see 
Eq.  (41)  in  Rlno  and  Fremouw,  1977). 


The  function  Q(q}  gives  the  shape  of  Che  SDF.  It  Is  normalized  so 
that  I q Q(q)  dq/(2Tt^)  - 1.  Thus,  for  a power-law  SDF  we  can  take 

2V-2 


Q(q) 


r(v  4-  i/2)/r(v 


1)<! 


, 2 , 2, 
Uq  + q ] 


V + 1/2 


(6) 


If  we  make  the  definition 


r(v  + l/2)/r(v  - 1)  (7) 

it  follows  that  for  K » q^.  The  parameter 

will  be  referred  to  as  the  strength  of  turbulence. 
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Now,  if  we  eubetltuto  Bq.  (5)  end  Eq.  (6)  into  Bq.  (2)  end  ebeuge 
verleblee,  the  reeult  ie 


Ng^(y)  ■ •«*=  9 G C 


Iq  Jo  <‘0 

“1 — T 

Cq:  + q 1 


(qy) 

m/2 


dq/2ic 


(8) 


where 


G - 


eh 


- B^/4  coe  8 

end  y le  replaced  by  f(dp^)  where 


(9) 


2 -* 


_* at  X X 


AC  - B /4 


<10) 


A discussion  of  Eq.  (8)  thet  describes  Its  relation  to  the  Briggs  and 
Perkin  (1963)  formulation  is  given  in  Rlno  and  Fremouw  (1977), 

The  integral  in  Eq.  (8)  can  be  evaluated  giving  the  result 

-(v-1/2)  v-i/2. 


850^7)  • “ec  6 G 


^-1/2 


2it  r(v  + 1/2) 


(U) 


where  K^(x)  is  the  modified  Bessel  function.  In  an  actual  experiment, 
we  TOuld  measure  a tenqwral  autocorrelation  of  the  form  where 

,-J/2 


'eff 


* 2 2 
Cv  - Bv  V + Av 

— J5 jeJz ja 


AC  - B /4 


and 


v^  - tan  e 


% 


(12) 


(13) 


^ Mb 

Finally,  v ■ relative  scan  velocity  at  the  ionospheric 

penetration  point  induced  by  both  the  source  motion  and  the  Irregularity 
drifts.  The  effective  scan  velocity  parameter  v^^j  cannot  be  larger  than 
v.  However,  Vg£f  can  be  substantially  smaller  than  v- for  exmnple,  if  the 
scan  direction  is  along  the  principal  irregularity. 
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ui'  1^  ieoti«ii£*ir  tb«  tmupoval  pomr  tpeetiruiB  of  i^imaai  wliich  is 
ikjatmA  1^  r.h«  «acpT««sioia 


«p(£)  - I *t>i  )I2idLiB5t)  A 4t 


m) 


«aai 


Sud^stitutlDis  Ef.  (11)  into  E*}.  <14 ) and  •vvluatlng  tlta  lutegsral  glvoa 


<Sp(£)  «*  sec  9 G C 

® *2/^  rOM-1/2)  V 


e££  Iq^  + (2«f }'* 


(15) 


glow,  recalling  the  ieflnltlon  of  the  pes^oaetc::*  [Eq.  (7)],  we  make 
the  analogpue  definition 


T 


2.2 

'e^ 


(L  sec  0)  G C 

9 


Viirtv) 

(2i<)^*  r<v+l/2) 


2V-1 


eff 


(16) 


~2v 

so  that  whsnsver  <2r£/v  ,.)  » q , 9(f ) » Tf  . As  a check  on  the 

eii  O' 

computations,  we  note  that  the  spectral  Index  of  the  one-dinienalonal 
phase  Spectrum,  .2Vt,  is  one  less  than  the  corresponding  spectral  Index  of 
the  three-dimens icrnal  Irregularity  spectrum,  2v+l  Isee  Eq.  (6)1, 

Equation  (16)  shows  that  the  phase  scf.ntlllatlon  level  depends  critically 
on  the  propagation  geometry,  particularly  through  G and  v ^r,  as  well  as 

»ia  ' 

the  relative  scan  velocity  v. 

To  continue,  we  note  that  a detrending  procedure  (cssentlelly  a 
high-pass  filter)  must  be  applied  to  separate  the  slow  TE€- induced  Crend- 
llke  phase  v-ariations  from  the  more  rapid  phase  scintillations.  It  Is 
this  detrending  procedure  that  dictates  the  smallest  measurable  temporal 
frequency  component  in  «p(f)-“8ay,  f^.  Since  we  have  found  no  systematic 
intrinsic  cutoff  cutoff  in  the  8pectr*.im,  we  must  conclvtde  that  the  in- 
equality <1^  always  holds.  It  follows  that  tUe  only 

unambiguous  characterization  of  the  phase  spectrum  Is  in  terms  of  the 
parameters  T and  p » 2v,  which  are  routinely  measured  in  the  Wideband  data 
reduction  (Fremouw  et  al.,  1978). 


The  consequences  of  this  fee;:  are  lo^rtent-,  Firsts  oo  absolute 
value  or  even  uppe:^'  bound  can  be  assigned  Lo  the  phase  variance.  The 
measured  phaie  variance  Is  reasonably  well  approximated  by  ^-he  foramla 


(6$b  ^ T 


d£ 


J j-pn  r(v.i/2)  _ 

® r(v) 


From  Eq.  (16)  and  Kq.  (17)  it  follows  that 


<5^^)  « X^(L  sec  6)  G • 


(18) 


Now,  if  we  take  the  limit  of  Eq.  (11)  as  y -•  0 or  Integrate  V(f)  over  ell 
frequencies,  we  obtain  the  ideal  phase  variance 


2,2 


<60^>  - i';X 


(L  sec  0)  G 


q^^^^  r(v-l/2) 
4«  r(v+l/2) 


(19) 


\dilch  la  the  conventional  rms  phase  expression. 


To  properly  Interpret  icncispheric  phase  scintillation  data,  Eq.  (18) 
must  be  used.  The  principal  difference  between  Eqs.  (18)  and  (19)  ir 
the  dependence  of  the  former  on  Indeed,  depends  critically 

on  altitude.  Hence,  the  measured  rms  phase  will  also  depend  on  altitude. 


To  sittomerlze,  an  unsmbiguous  characterization  of  phase  scintillation 
data  can  be  obtained  only  in  terms  of  the  spectral  strength  parameter  T 
and  the  spectral  Index  p.  If  we  con  estimate  the  anisotropy  and  drift 
of  the  Ionospheric  irregularities,  Eq.  (16)  can  be  used  to  estimate  the 
strength  of  turbulence  \dilch,  as  with  T Itself,  is  presently  th-i  only 
unambiguous  average  parameter  that  con  be  used  to  characterize  the 
spectral  strength  of  the  ionospheric  irregularities. 


We  have  assumed  In  this  analysis  that  diffraction  effects  in  the 
phase  data  are  negligible.  In  Fremouw  et  al.  (1978),  it  is  slkown  that 
the  measured  rms  phese,  when  corrected  for  the  finite  reference  fre- 
quency [Eq.  (11)1,  varies  linearly  with  wavelength.  We  shall  now  show 
that  under  conditions  of  weak  scatter  the  detailed  structure  of  the  phase 
Itself  scales  with  frequency  as  predicted  by  Eq.  (1)  to  periods  shorter 
than  one  second. 
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In  Figure  i ve  Bh<y»  a typical  V8S  phase  scintillation  record.  The 
raw  phase  data  have  been  detrended  to  resxive  phase  variations  with  periods 
greater  than  10  s (f^,  - 0.1  Ha).  Beiore  0920  VT  the  8^  scintillation 
index  is  leas  than  0.A,  In  Figure  2 we  show  on  an  espanded  scale  the 
differertces  betweaa  the  phase  at  the  J.ndlceted  lower  frequencies  and  the 
scaled  phase  at  the  Indicated  higher  frequencies.  The  designations  UL3 
and  UV3  denote,  respectively,  the  lowest  and  hipest  of  seven  equiapaced 
UHF  frequencies.  One  can  see  that  prior  to  0920  UT  the  large-scale  phase 
variations  are  completely  suppressed,  and  therefore  unaffected  by  dif- 
fraction. 


FIGUne  1 SEGMENT  OF  DETRGNDED  PHASE  DATA  (fc  • 0.1  Hr)  FROM 
•VIDEBAND  PASS  6-33  RECORDED  AT  POKER  FLAT,  ALASKA 

To  demonstrate  this  quantitatively,  in  Figure  3 we  show  the  VHF 
phase  spectrum  and  the  spectrum  of  the  VHF-UU3  phase  difference.  Below 
1 Hz  there  is  essentially  20  dB  of  cancellation.  Beyond  30  Hz,  the  phase 
difference  spectrum  is  ~ 3 dB  above  the  VHF  phase  spectrum  Indicating 
total  decorrelation.  We  see  that  the  diffraction  effects  are  largely 
confined  to  frequencies  greater  than  1 Hz,  However,  the  dominant 
spectral  content  comes  from  frequencies  below  1 Hz,  which  explains  why 
the  ms  phase  follows  the  linear  wavelength  dependence  so  accurately 
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FIGURE  2 DIFFERENCE  BETWEEN  MEASURED  PHASE  AT  LOW  FREQUENCY  AND  SCALED 
PHASE  FROM  HIGHER  FREQUENCY  FOR  SEGMENT  OF  POKER  FLAT  PASS  6-36 


0.1 


1.0 


10 


FREQUENCY  — Hi 
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FIGURE  3 SPECTRAL  DENSITY  OF  VHF  PH/JE  AND  VHF-UU3  PHASE  DIFFERENCE 
FOR  20-s  DATA  RECORD  FROM  POKER  FLAT  PASS  6>36  AT  ABOUT 
0919  UT 


We  also  conclude  from  Figure  3 that  errors  In  T due  to  diffraction, 

idilch  act  to  flatten  the  spectrum,  are  less  severe  than  the  corresponding 

errors  In  p.  Strong  scatter  effects,  Mhlch  are  evident  In  Figure  2(a) 

after  0920  UT,  also  generate  mainly  high-frequency  components.  Ultimately, 

-2 

strong- scatter  effects  will  drive  the  phase  spectrum  toward  f . By 
careful  processing  and  data  Interpretation,  we  can  usually  rocognlse  and 
avoid  erroneous  T and  p values. 
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•Jt: 


Ill  INTBMSm  SCIKTIUATIOH 


In  this  section  «•  shell  develop  the  form  of  the  8^  scintllletlon 
Index  In  the  limit  of  en  Infinitely  lerge  outer* scale  cutoff.  'j1:1s 
limit  Is  well  defined  because  Tresnel  filtering  suppresses  the  contri- 
bution of  large-scale  coiq>o<^Anta.  The  results  take  a fairly  slaqtle 
analytic  form  that  properly  accounts  for  the  propagation- angle  dependence 
of  the  scintillation.  In  Sufenach'a  (1975)  formulas,  which  retain  an 
explicit  dependence  on  the  outer  scale,  the  angle  dependence  is  Introduced 
In  an  ad  hoc  manner. 

We  shall  characterise  the  amplitude  scintillation  by  the  second- 
order  Dements  of  intensity.  However,  uMer  the  asstansd  weak-scatter 
conditions,  the  atoplltude  and  log-amplitude  are  slnply  related.  For 
example,  the  scintillation  Index  S^,  which  Is  the  normalized  rms  intensity, 
is  twice  the  rms  amplitude  index.  The  Intensity  autocorrelation  function 
corresponding  to  Eq.  (2)  Is 


R_(A3  ) ■ L ® 

X 8 6 


X sin^  (h(lC)z)  cos  (tC  . /!lsp  ) 


dlt 
(2«)^ 


•C) 


(20) 


where 


h(l?)  - . ten^  e( 


(21) 


and 


Xz  sec  6 
4Tt 


(22) 


The  square  of  the  scintllletlon  Index  Is  derived  from  Eq.  (20) 

by  lotting  £4}  * 0.  If  ve  substitute  Kqs,  (5)  and  (6)  Into  Eq.  (20), 

8 

use  Eq.  (7),  and  then  change  variables,  we  obtain  the  simpler  expression 
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2 2 2 

jf  " 4r  A (L  8«c 
4 e 


JJ  [qj  + (aV  + 


+ C'qp] 


v+1/2 


X sln^  t(q^  + qy)Zl 


aq 

2ff 


2jt 


(23) 


Th®  coefficient*  A^,  end  ere  related  to  A,  B,  and  C by  Eq.  (51)  In 
In  Rlno  and  Fraeiouw  (1977)  where  It  1*  shown  that  Eq.  (23)  Is  equivalent 
to  the  Brlggs-Parkln  formula  as  genarallxad  by  Singleton  (1973)  [see  also 
Rufenach  (1975)  and  Costa  and  Ktlley  (1976)]. 


As  It  stands,  Eq.  (23)  cannot  be  evaluated  analytically.  Rufenach 
(1975)  and  Costa  and  Kelley  (1976)  have  used  a sllg^htly  modified  form  of 
Eq.  (23),  which  can  be  evaluated  when  v - 1.5.  However,  from  the  discus- 
sion In  Section  I,  it  follows  that  for  Ionospheric  sclntllletlon, 

2 

« 1.  Thus,  It  la  a good  approximation  to  take 


2 2e 

S.  • 4r  (L  sec  0)  C 
^ © s 


V i/-ff  ^ 

(A'q^  + ^ C^qJ) 


(24) 


which  Is  the  limiting  form  of  Eq.  (23)  as  q^  -•  0.  Hereafter,  we  shall 
Jenote  the  double  Integrrl  in  Eq.  (24)  by  I, 


To  evaluate  I.  we  fl>:8t  perform  a rotation  of  coordinates  to  remove 


the  q q term.  The  result  xt 
X y 


■I 


ab  sin^  (q^)  dq  dq  /(2;t)^ 
X y 


(A  Rjj  + C qy) 


V + 1/2 


(25) 


where 


and 


ifA'  + C'  + d"] 

(26a) 

\lk'  + c'  + d'] 

(26b) 

/-.  / in 'a  2 #2 

V(A  - C ) + B 

(26c) 

We  note  that  A*  ^ C^, 
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T"a™ 


4 ■ ’ 

- . * * , ' f - ■ 


Mov,  1 is  s«f»srs1»U  such  ttlist:  sfttr  s strltis  pi  vmfilpplmHemPf 
I - f •In*  «i’)  *t  f • •'■•V  ( 

' * Jl  U'  - <*'-0l  •in®  *1 

Ihs  Intsgrsl  over  q is  wsXl  koown  for  ccsttsritkg  by  Isotropic  itrigo* 
larltisa*  The  intsgrsl  ovst  4 properly  accounts  for  the  gsometricsl 
offsets  of  anisotropic  irrsgularitlsft.  In  contrast,  Sufsnach  (l$7j) 
introduesd  on  sd  hoc  oultlplicsti^  gsometrlcal  factor  to  account  for 
anisotropic  irregular Itiss*  By  evaluating  the  intsgrsl  over  q and 
substituting  Bq.  (27)  into  Eq,  (24)  ve  have 

- .2.2 „ .V-I/I  [•r(^)«.  U(I-2v)/4l] 


a.  - 4r^X  (L  ssc  0)  CZ 

M‘  G 8 


2r  2<5  - 


shsra 


i^(v+l/2) 


it/2 

lb  I d0 

•1/2)  I r r'  9 nv  + 172 

jo  1 - A...:. g.  Bln^  ^ 

I j 


The  integral  in  Bq.  (29)  can  be  evaluated  in  terns  of  the  hyper- 
gsonetrical  function  GradBhte}^  and  Ryshik,  1965; 
Formula  3.681).  To  avoid  convergence  problems  when  A*'  » c',  however, 
we  have  applied  the  transformation 

2Fj(q',P;v;*)  « (l-z)''"“'^2^l(''-®»  ''■P»  ''i  • (' 
Making  the  appropriate  substitutions,  7 can  be  evaluated  as 


vS’-c-''  ® 


-j(l/2  - V,  1/2,  1, 


Now,  r’C'  + 1/2)1.  By  direct 

computation  from  Bq.  (27)  when  a » 1,  it  can  be  shown  that  this  limit 


is  correct. 


To  lunMrlao,  Si|s.  (26),  (28),  and  (31)  can  bo  uoad  to  avaluato  tha 
scintillation  indax  undar  conditions  of  waak  scattering  In  tanas  of  0^, 
V,  z,  and  tha  propagation  gaoiaatry.  »a  note  that  tha  wavelength  depen- 
dence of  S.  Implied  by  Eq.  (28)  Is  S.  « xhus,  if  v - 1.5, 

^ - 3 1 5 

which  corresponds  to  <P(f)  f , « X * , the  nominal  wavelength 

dependence  typically  reported  for  Ionospheric  data. 

By  using  Eq.  (16),  we  can  write  Eq.  (28)  In  the  equivalent  form 

C(v)  (?/«)  (32 

where 


C(v) 


cos 


(rt<l-2v)/4]  r(v+i/2)  (2rt) 
(5  - 2v)72 


2v 


r(v)  2 


(33) 


If  wc  note  that  the  units  of  T are  radians  squared  per  hertz  raised  to 

the  2V  + 1 power,  It  is  easily  verified  that  Eq.  (32)  Is  dimensionslly 

consistent.  We  also  see  from  Eq.  (32)  that  for  a fixed  rms  phase  level, 

varies  inversely  with  veff,  which  increases  with  increasing  height. 

Thus,  while  the  factor  to  Increase  S.  with  increasing  height, 

-(2v-l) 

the  decreasing  factor  dominates  and  actually  decreases. 

As  we  shall  see,  this  effect  Is  important  In  interpreting  the  Wideband 
satellite  data. 
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IV  APPLICATION  TO  WIDEBAND  SATELLITE  DATA 


A.  General 

We  shall  now  apply  the  phase-screen  model  to  the  interpretation  of 
Wideband  satellite  data.  As  discussed  In  Fremouw  et  al.  (1978),  we 
routinely  measure  the  phase  SDF  at  VHF  (137  MHs)  and  UHF  (378.6  MHs). 

The  spectral  estimates  are  smoothed,  after  which  a log- linear  least- 
squares  fit  Is  applied  to  determine  the  spectral  strength  parameter  T and 
the  spectral  Index  p.  Thus,  ^^(f)  ^ Tf  ^ over  the  significant  portion  of 
the  phase  SDF.  The  frequency  limits  for  the  fit  are  carefully  chosen  to 
minimise  noise  contamination  and  detrend  filter  effects. 

We  shall  first  apply  Eq.  (17)  to  compare  the  calculated  and  measured 
rms  phase.  If  the  two  results  agree,  we  can  be  confident  that  the  phase 
SDF  Is  Indeed  accurately  modeled  by  the  power- law  form.  Thus,  this  first 
seep  is  mainly  a consistency  check  for  the  basic  parameters  T and  p. 

The  next  step  is  to  apply  Eq.  (32)  to  compute  S^  for  comparison  with 
its  measured  value.  To  evaluate  Eq.  (32),  however,  we  must  specify  the 
height  of  the  phase  screen,  z,  the  relative  scan  velocity,  v,  and  the 
anisotropy  parameters  a,  b,  6,  ^ere  a Is  the  euclal  ratio  along  the 
magnetic  field,  b is  the  axial  ratio  transverse  to  the  magnetic  field, 
and  fi  is  the  orientation  of  the  transverse  axis  such  chat  S 0 for 
geomagnetic  east-west  sheets  (see  Rlno  and  Fremouw,  1977). 

In  our  routine  summary  analysis  we  calculate  the  satellite  component 
of  the  relative  scan  velocity  and  the  propagation  angles  at  two  reference 
altltudes--namely,  110  km  and  350  km.  Thus,  if  the  satellite  component 
doml.nates  the  Irregularity  drifts  we  can  compute  at  E-  and  F-reglon 
altitudes  for  different  anisotropy  parameters.  When  we  achieve  a good 
fit  to  the  data,  we  can  remove  the  geometrical  factors  in  Eq.  (16)  to 

estimate  the  strength  of  turbulence  times  the  layer  thickness  LC  . 
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Becruce  we  have  no  direct  mfians  of  determining  the  effective  layer 

thickness,  LC  Is  the  most  basic  measure  of  the  Irregularity  strength. 

For  display  purposes,  however,  we  shall  specify  a representative  layer 

thickness  from  which  an  actual  value  for  C can  be  obtained.  The  inter- 

8 

pretation  of  will  be  discussed  in  Section  V. 

B.  Equatorial 

Before  considering  Individual  passes.  It  Is  useful  to  look  at  the 
average  behavior  of  the  spectral  Index  p.  As  we  noted  In  Section  II, 
the  measured  value  of  p Is  sensitive  to  noise  contsuiluatlon  under 
conditions  of  weak  scattering  and  to  diffraction  effects  under  conditions 
of  strong  scattering.  To  demonstrate  these  effects,  we  have  plotted 
the  average  value  of  p versus  S^  for  a representative  set  of  passes. 

The  Ancon  data  are  shown  in  Figure  4.  The  measured  p Index  achieves 
a maximum  value  near  but  slightly  less  than  3 for  values  between  0,4 
and  0.6.  Noise  contamination  end  diffraction  effects  cause  the  respective 
reductions  of  the  measured  p values  frr  weak  (S^  < 0,4.  and  strong 
(S^  > 0.6)  fading  levels.  The  wavelength  dependence  of  under  conditions 
of  weak  scattering  provides  an  independent  means  of  verifying  the  value 
of  p. 


FIGURE  4 AVERAGE  VALUE  OF  p FOR  S4  VALUES  WITHIN  THE 

INTERVALS  0.2  n < S4  < 0.2  (n  * 1)  AS  DERIVED  FROM 
ANCON  DATA.  The  UHF  and  VHF  curves  are  computed 
separately. 
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In  Figure  5 we  show  a scatter  diagram  of  Ancon  values  measured 
at  VHF  and  UHF.  The  paucity  of  data  between  3^  * 0.2  and  ■ 0,5  Is 
due  to  the  tendency  of  the  equatorial  VHF  scintillation  to  be  either 
weak  or  strong  (Llvlr^stony  1978).  In  any  case,  the  curve,  ^Ich 

corresponds  to  p » 3,  fits  the  data  quite  well.  The  fit  Is  more  striking 
In  Figure  6,  where  we  have  plotted  at  UHF  against  at  L band. 


FIGURE  5 SCATTER  DIAGRAM  OF  S4  AT  VHF  vs  S4  AT  UHF  FOR 

A SUBSET  OF  THE  ANCON  PASSES  USED  IN  GENERATING 
THE  CURVES  IN  FIGURE  4 


The  average  p values  for  the  Kwajaleln  data  show  behavior  similar  to 
to  that  of  the  Ancon  data  (see  Figure  7).  However,  the  maximum  average 
p value  falls  distinctly  below  p = 3.  The  corresponding  shallower  wave- 
length dependence  of  is  consistent  with  a smaller  p Index  as  shown  In 
Figures  8 and  9 where  we  have  plotted  the  scatter  diagrams  for  VHF 
versus  UHF,  and  UHF  versus  L band.  Thus,  It  appears  that  the  phase  spectra 
obtained  from  the  Kwajaleln  data  are  syetematlcally  flatter  than  the  phase 
spectra  obtained  from  the  Ancon  data. 
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FIGURE  6 SCATTER  DIAGRAM  OF  S4  AT  UHF  vt  S4  AT  L-BAND 
FOR  A SUBSET  OF  THE  AMCOH  PASSES  USED  IN 
GENERATING  THE  CURVES  IN  FIGURE  4 


In  both  the  Ancon  data  (Figures  5 and  6)  and  the  Kwajaleln  data 
(Figures  8 and  9)  saturation  of  the  Intensity  scintillation  apparently 
occurs  for  an  value  less  than  unity,  (/n  value  of  unity  corres- 
ponds to  Rayleigh  fading.)  This  effect  Is  due  to  our  detrending  operation, 
which  Is  applied  separately  to  Intensity  and  phase.  The  strong  fading 
data  from  the  equatorial  stations  evidently  contained  significant  Fourier 
components  beyond  the  10  s cutoff  of  the  detrend  filter.  This  will  not 
be  a problem  for  our  analysis  here  because  ve  are  considering  only  weak- 
scatter  data  In  this  report. 

In  Figure  10  we  show  the  phase  and  Intensity  scintillation  data  for 
a disturbed  nighttime  Ancon  Wideband  pass  together  with  a set  of  theoret- 
ical calculations  of  and  Consider  first  the  rms  phase.  We  see 

that  V « 1.3  gives  as  good  a fit  to  the  data  as  the  v * 1,5  curve,  which 

.3 

corresponds  to  an  f phase  SDF.  As  we  have  already  noted,  diffraction 
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FIGURE  7 AVERAGE  VALUES  OF  p FOR  S4  VALUES  WITHIN  THE  INTERVALS 
0.2  n < S4  < 0.2(n  + 1)  AS  DERIVED  FROM  KWAJALEfN  DATA 

effects  cause  a flattening  of  the  phase  SDF,  which  explains  the  dis- 
crepancy between  the  measured  and  calculated  nns  phase  after  ~ 0345  UT. 

The  F- region  curve  In  Figure  10  for  v ■>=  1.3  comes  cloeer  to  the 

measured  curve  than  the  corresponding  v « I.5  curve.  However,  in 
the  regime  where  S^  Is  small  and  we  expect  the  weak-scatter  theory  to  be 
applicable,  the  theoretical  curve  gives  a result  that  Is  consistently  too 
large.  As  we  noted  at  the  end  of  Section  III,  however,  raising  the 
equivalent  phase  screen  helghc  lowers  for  a specified  rms  phase  level. 

Thus,  we  believe  that  the  discrepancy  In  Flpire  10  la  due  to 
the  350-km  reference  altitude  being  too  low.  Unfortunately,  the  pr^a- 
gatlon  angles  are  only  cemputed  for  two  heights.  Nonetheless,  In 


Figure  10  we  tihow  the  cexculetions  for  the  E- region  reference  altitude. 
Again,  the  large  eatlnatea  for  the  E-reglon  ere  due  to  the  smaller  scan 
velocities  [see  Eq.  (32)].  Finally,  the  large  overestimates  of  when 
Is  greater  than  ~ 0.4  are  due  to  the  fact  that  the  weak-scatter  theory 
does  not  allow  for  multiple  scattering. 

To  Illustrate  the  axlal-ratlo  dependence  of  we  have  used  the 
geometry  for  Ancon  Pass  29*7,  but  have  assigned  a fixed  value  to  C . 

Tne  corresponding  values  for  Isotropic  Irregularities  and  axial 
ratios  of  5:1  and  100:1  are  shown  In  Figure  11.  Both  the  variation  of 
with  changing  propagation  geraietry  and  the  magnitude  of  decrease 
with  Increasing  axial  ratio  until  the  axial  ratio  exceeds  10:1.  Beyond 
10:1,  the  index  for  the  equatorial  geometry  does  not  eidiibit  an  axial 
ratio  dependence.  Ttte  diffraction  is  then  effectively  two-dimensional, 
as  discussed  In  Section  III. 

To  complete  the  equatorial  examples  from  Ancon,  we  have  selected 
three  additional  passes  and  applied  the  nominal  F*reglon  geometry  with 
V « 1.4  and  a 100:1  axial  ratio.  The  results  are  shown  In  Figures  12, 

13,  and  14.  In  Figure  12,  the  calculated  values  fall  slightly  below 
the  measured  values.  For  such  low  values,  noise  contamination  Is  a 
possible  explanaticn  (see  Appendix).  However,  Increasing  the  v value 
and/or  lowering  the  reference  altitude  will  tend  to  Increase  the  model 
values  as  we  have  noted. 

In  Figure  13,  the  measured  values  are  significantly  larger  th«m 
those  shown  in  Flgu.'s  12.  Thus,  we  expect  the  model  calculations  to  over» 
estimate  S^.  There  Is  a gener.  1 tendency  for  the  measured  phase  to  fall 
below  the  theoretical  curve,  which  Is  evidently  the  diffraction  effect 
we  have  already  noted.  Nonetheless,  the  C values  should  be  roughly 
correct,  and  Indicative  of  the  turbulence  levels  required  to  produce 
significant  gigahertz  scintillation.  The  data  set  shown  In  Figure  14 
falls  between  the  extremes  shown  In  Figure  12  and  13.  The  overall  fit 
In  Figure  14  Is  quite  good. 

Turning  now  to  the  Rwajaleln  data.  In  Figure  15  we  show  a typical 
disturbed  nighttime  pass.  The  rms  phase  calculations  for  v • 1.25  clearly 
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fit  the  Tins  phase  data  better  than  the  calculations  using  v ■ 1.5.  As  with 
the  Ancon  data,  however,  there  is  a tendency  for  the  F«region  curve  to 
slightly  overestimate  the  data  which  we  have  attributed  to  increased 
layer  height.  Thus,  the  only  systematic  difference  between  the  Ancon 
and  Kvajaleln  data  that  can  be  ascertained  from  the  first-order  moments 
is  the  lower  phase  spectral  index— p ~ 2.5  for  Kwajaleln  versus  p 2.8 
for  Ancon.  This  difference  implies  somewhat  sharper  phase  gradients  in 
the  Kwajaleln  data  than  in  the  Ancon  data. 

As  a final  conanent  on  the  equator  data,  we  have  not  made  direct 
computations  at  Irband  because  the  I^band  phasu  data  are  not  routinely 
Fourier  analysed.  On  the  other  hand,  we  have  verified  the  wavelength 
dependence  of  between  UHF  and  L-band  in  Figures  6 and  9.  Moreover, 
even  if  UHF  is  near  saturation,  the  phase  scintillation  ard/or  T can  be 
frequency  scaled  with  u high  degree  of  accuracy.  Thus,  the  self- 
consistent  Intensity  and  phase  calculations  at  UHF  can  be  accurately 
extrapolated  to  L-band.  An  example  is  presented  in  Section  V. 

C.  Auroral 

The  interpretation  of  the  auroral- zone  data  Is  complicated  by  the 
rapidly  changing  propagation  geometry.  Moreover,  auroral-zone  scintil- 
lations are  typically  associated  with  extreme  variations  in  perturbation 
strength,  albeit  at  substantially  smaller  levels  than  the  equatorial 
data.  As  with  the  equator  data,  we  begin  hy  considering  the  general 
behavior  of  the  spectral  Index. 

A plot  of  p versus  for  a representative  sampling  of  Poker  Flat 
passes  Is  shox^n  In  Figure  16.  Here  the  average  spectral  index  obtained 
from  the  VHF  phase  spectra  is  less  than  2.5.  A set  of  disturbed  passes 
were  used  to  obtain  the  UHF  values,  which  tend  to  give  a slightly  larger 
p Index.  Tnus,  there  may  be  a systematic  steepening  of  the  auroral  phase 
spectra  with  increasing  perturbation  strength,  or  associated  with  those 
events  that  produce  significant  UHF  scintillation. 

In  general,  it  appears  that  the  auroral  data  tend  to  show  even 
steeper  phase  gradients  than  the  equatorial  data  from  Kwajaleln. 
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FIGURE  16  AVERAGE  VALUE  OF  p vs  S4  AS  IN  FIGURES  4 
AND  7.  The  VHF  data  are  derived  from  SO  repre- 
sentative Poker  Flat  passes.  The  UHF  data  are 
derived  from  four  disturbed  Poker  Flat  pestM. 

However,  the  auroral  data  heve  shown  considerably  more  variability  In 
the  measured  p Index  and  the  overall  structure  of  the  scintillation. 

This  is  evident  in  the  scatter  diagram  of  S4  values  at  UHF  versus  S4 
values  at  VHF  shown  in  Figure  17.  The  curve  corresponds  to  a p 

index  value  of  2,  which  is  more  typical  of  reported  in- situ  measurements 
than  phase  scintillation  measurements. 

As  a first  detailed  example  of  auroral- zone  data,  we  have  selected 
a nearly  overhead  pass  with  an  isolated  scintillation  enhancement  near 
the  point  or  closest  approach  to  the  magnetic  zenith.  Such  events  are 
coimaonly  observed  under  conditions  of  moderate  auroral  activity.  The 
pattern  of  these  events  in  passes  to  the  east  and  west  of  the  station  is 
such  that  they  tend  to  occur  where  the  propagation  path  lies  within  the 
local  L-shell.  Thus,  we  have  hypothesized  a sheet-like  anisotropy  for 
the  irregularities. 

In  Figure  18  we  show  the  data  and  theoretical  calculations,  for  east- 
west-aligned  sheets  with  a b = 10  (10:10:1  notatlonally)  at  E-  and  F- 
reglon  reference  altitudes.  The  rms  phase  calculation  using  v « 1.25 
fits  the  data  very  well  except  near  the  localized  enhancement.  Here  the 


FIGURE  17  SCATTER  DIAGRAM  OF  S4  MEASURED  AT  VHF  vc 
514  MEASURED  AT  UHF  FOR  THE  FOUR  DISTURBED 
PASSES  USED  IN  GENERATING  THE  UHF  CURVE  IN 
FIGURE  1(1 

phase  SDF  steepens  such  that  p » 3.  We  belie\'e  this  effect  Is  due  to  the 
nonstatlonarlty  Induced  by  the  rapidly  changing  propagation  geometry. 
Indeed,  a purely  geometrical  effect  should  produce  no  change  In  the 
spectral  shape. 

The  E-reglon  and  F- region  S^  curve-  for  10:10:1  sheets  tend  to 
bracket  the  data  near  the  enhancement.  Hence,  If  the  10:10:1  sheet  model 
Is  correct,  the  equivalent  phase  screen  must  be  placed  between  the  E> 
and  F-reglon  reference  altitudes.  We  note,  however,  that  both  calculations 
underestimate  the  measured  S^  values  at  the  extremes  of  the  pass.  Rod- 
llke  irregularities,  however,  fit  the  data  quite  wellc 
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We  have  hypothesized  that  the  sheet- like  structures,  which  are 
Independently  verified  by  analysis  of  the  Wideband  spaced-recelver  data 
(Rlno  and  Livingston,  1978),  are  confined  to  the  region  of  the  diffuse 
(continuous)  and  discrete  aurora.  Outside  these  precipitation  regions, 
the  Irregularity  structures  are  more  nearly  rod-like,  'fo  Illustrate  this 
posslbili'.ty,  we  also  show  the  theoretical  curve  for  8:4:1  sheets  In 
Figure  18.  The  8:4:1  curve  underestimates  the  enhancement,  but  fits  the 
extremes  quite  well. 

Now,  It  has  been  noted  from  the  earliest  Wideband  data  (Fremouw  et 
<il,,  J978''  that  the  auroral-zone  data  show  a large  number  of  events  In 
which  the  phase  scintillation  level  Is  much  larger  than  the  corresponding 
Intensity  scintillation  level,  even  though  the.  latter  remains  In  the 
weak-scatter  regime.  The  geometrical  enhancement  produces  such  an  effect. 
Tlie  rms- phase- to- ratio  at  the  rms  phase  peak  for  the  data  In  Figure  18 
approaches  4.0.  We  shall  see,  however,  that  there  are  other  events  that 
produce  large  rms  phase  enhancements  without  a proportionate  enhancement 
In  S^. 

In  Figure  19  we  show  an  exacxple  of  a disturbed  low-elevatlon  pass. 

The  rms  phase  data  fit  the  v = 1.25  curve  better  than  the  v » 1,5  curve. 
Near  0912  UT,  diffraction  effects  tend  to  reduce  the  measured  rms  phase 
below  the  predicted  value  Just  as  in  the  equator  data.  The  predicted 
values  for  the  F-regio.  10:10:1  sheet  model  fit  the  measured  data  very 
well  in  the  region  between  0909  UT  and  0911  UT.  Before  0909,  however, 
the  theoretical  F-region  calculations  underestimate  the  measured  values 
of  and  the  E-reglon  calculations  give  a better  fit.  For  completeness, 
we  also  show  the  v = 1,5  F-region  curve,  which  clearly  gives  a poorer  fit 
to  the  measured  curve  as  well  as  to  the  rms  phase  data. 

The  change  In  structure  that  occurs  at  0909  UT  is  evidently  due  to 
an  Increase  In  height  of  the  corresponding  Irregularities.  Recall  that 
for  a given  phase  scintillation  level,  the  self-consistent  value 
decreases  with  increasing  height,  and  vice  versa.  This  implies  that 
increasing  the  layer  height  tends  to  enhance  the  ratio  of  phase  to 
amplitude  scintillation.  Unlike  the  geometrical  enhancement,  however,  a 
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change  in  Irregularity  hei^t  would  produce  a roughly  equivalent  increate 
in  the  equatorial  and  auroral  date.  We  noted  one  such  event  in  Figure  10. 

In  Figure  20  we  show  another  example  of  a elevation  paaa,  but 

with  a strong  perturbation  present.  With  the  exception  of  the  regions 
of  large  values^  the  F-region  calculations  with  v * 1.25  and  v > 1.5 
fit  the  data  reasonably  veil.  However,  between  1052  and  1054  UI  the 
V <*  1.25  curve  underestimates  the  xms  phase  while  the  v > 1.5  curve  ever- 
estimates  the  rms  phase.  Allowing  for  diffraction  effects,  therti  is 
evidence  here  of  a systematic  change  in  the  phase  spectrum  within  the 
pass^ 

In  Figure  21  we  show  a more  severe  enample  of  such  an  effect.  We 
see  that  the  rms  phase  enhancements  are  badly  underestimated  by  the 
theoretical  calculations,  particularly  between  0948  UT  aud  0952  UT.  We 
believe  thet  this  is  due  to  rapidly  changing  perturbation  structures 
that  do  not  readily  admit  a homogeneous  spectral  representation,  rather 
than  to  some  fundamental  change  in  irregularity  developnient  leading  a 
non- power-?. aw  specti'al  distribution.  Either  way,  the  effect  is  to  produce 
large  phase  excursions  without  a proportionate  increase  in  S^. 

Such  events  are  not  uncommon  in  the  auroral  data.  However,  while 
they  are  not  properly  characterised  using  signal  moments  calculated 

under  the  weak-scatter  theory,  the  measured  anisotropy  and  relative 

( 

pattern  drifts  obtained  from  the  spaced  receiver  data  arc  veil  behaved 
(see  Rlno  and  livingston,  1978).  It  is  from  such  an  anal3r8lR  that  the 
E-region  model  with  a - 8 and  b = 4 ir  Figure  21  was  deduced.  It  is 
interesting  that  in  spite  of  the  poor  fit  to  the  rms  phase  data,  the 
E- region  mdcl  calculations  fit  the  data  reasonably  well. 

To  sumnarlre  the  auroral  data,  we  have  sliown  that  as  long  as  the 
rms  phase  data  fit  the  power- law  model  and  is  less  than  0.4,  the 
weak-scatter  pov^er-iaw  phase-screen  model  gives  self-consistent  results 
for  rms  phase  and  intensity  with  an  appropriate  choice  of  irregularity 
(phase- screen)  height  and  irregularity  anisotropy.  Indeed,  the  model  is 
particularly  sensitive  to  irregularity  height  and  detailed  anisotropy. 
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FIGURE  20  VHF  DATA  FROM  POKER  FLAT  PASS  6-15  SHOWING  EVIDENCE  OF 
CHANGE  IN  PHASE  SPECTRA 


We  have  taken  F-reglon  irragularitlaa  with  a 10:10tl  ehaat»llke  | 

anisotropy  as  representative  of  auroral-sone  Irregularities.  Howevery 
we  have  not  attempted  to  demonstrate  unaqulveeally  here  that  the  sheet 
model,  as  opposed  to  a slnple  rod- like  anl8otrop\ , is  correct.  A 
separate  report  will  address  the  Irregularity  anisotropy  lu  detail.  A | 

point  to  be  made  here,  however.  Is  that  Wherever  a localieed  geoomtrlcal  | 

enhaacement  occurs,  the  ms  phase  Increases  more  rapidly  than  the  increase  j 

in  S^. 

We  have  also  shown  non- geometry- associated  events  that  produce 
large  rms  phase  enhancements  without  a proportionate  increase.  Such 
events  show  pronounced  departures  from  the  simple  power- law  model  we  i 

have  eiiq>loyed  in  our  analysis.  Hence,  the  theory  is  inappropriate  in  | 

such  cases.  We  beli'^ve  that  such  events  are  representative  of  undeveloped  ^ 

turbulence.  A smooth  ai  roral  arc  is  representative  of  a structure  that  | 

can  produce  large  localized  phase  enhancements  without  a proportionate  | 

Increase  in  S..  It  Is  clear  that  such  structures  are  not  appropriately  | 

^ Si 

characterized  by  second-order  moments.  | 

1 
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V PISCUS8I0» 

In  this  raport  va  have  first  davalopad  rma  phase  and  fomulas 
based  on  a power* law- phase- screen  nodal  with  an  arbitrarily  large  out%r- 
scale  cutoff r The  fonmiles  allow  fully  for  the  angle  dependence  of  the 
signal  moments  In  a highly  anisotropic  medium.  We  have  shown  that  the 
power- law- phase- screen  model  gives  an  accurate  self-consistent  description 
of  the  phase  and  Intensity  scintillation  under  condltlont.  of  weak 
(S^  « 0,4)  intensity  scintillation. 

The  height  of  the  equivalent  phase  screen  and  the  anisotropy  of  the 
irregularity  are  the  only  free  parameters  In  the  theory.  We  have  found 
that  the  parameter  sensitivity  of  the  model  Is  such  that  re.-^sonable 
bounds  can  be  put  on  the  height  and  anisotropy.  We  1 ave  not,  however, 
tried  to  carefully  determine  the  morphology  of  the  Irregularity  height 
here  oecause  this  can  be  better  accomplished  by  uslug  the  Interferometer 
data.  That  analysis  is  being  reported  separately  (Rlno  and  Livingston, 
1978). 

Our  main  purpose  in  presenting  the  analysis  in  this  report  is  to 
demonstrate  the  adequacy  of  the  phase-screen  model  for  translating  the 
routinely  measured  Wideband  summary  yarameters  to  invariant  irregularity 
strength  measures.  Such  transli.tlons  arc  necessary,  for  exmple,  to 
predict  phase  perturbation  levels  that  might  adversely  affect  advanced 
survelll  iXV'fi  satellite  systems.  They  also  enable  us  to  predict  intensity 
scintillation,  particularly  the  equatorial  gigahertz  scintillation  that 
can  adversely  affect  satellite  connunicatlon  system... 

To  demonstrate  how  the  model  can  be  used,  let  us  consider  a phase 
screen  at  350  km  with  a signal  source  at  three  times  this  height  directly 
overhead  at  the  geomagnetic  equator.  We  can  then  use  Eq.  (27)  to  predict 
as  a function  of  C^.  Since  we  have  used  200  km  for  the  layer  thickness 
in  deriving  rapresentaflve  values  from  our  data,  we  shall  also  apply 
that  value  here.  For  lllustretlve  purposes,  we  have  used  isotropic 
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irr«gul«rltlei  and  100:1  rods.  JJhu®,  An  Flfore  22 (n)  we  hav«  olotted 
curve#  at  1 GH*  \^r#u8  for  v - 1,25  and  the  corresp ->tuilng  curve® 

In  Figure  22(b)  for  v ■ 1.5,  ua  aae  that  levels  approaching  10^^ 
correspond  to  significant  levels  of  glgahert®  scintillation. 

Our  entire  analysis  is  based  on  thfi  assumption  that  tha  three* 

dimensional  irregularity  SDF  has  the  form  C^/lq^  + The 

d o 

corresponding  spatial  correlation  function  has  the  form 
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It  is  easily  ahown  that  (ZiN  ) 
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is  consistent  with  Eq,  (7). 


as  determined  by  evaluating 
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Now,  if  an  In- situ  probe  scans  the  medium,  the  spectral  density 
function  that  characterizes  the  probe  output  is  obtained  by  Fourier 
transforming  Eq.  (34)  with  y replaced  by  v^6t,  where  v^  is  the  effective 
probe  velocity  analogous  to  Eq.  (12).  The  resulting  SDF  i® 


Cg  r(v*l/2)  ^ 

4jtV(v+l/2)  Vj[q^  + (2rtf/Vj)^]''"^^^ 


(35) 


When  2«f/vj  » q^,  tpp(f)  * T^f^''"^.  Thus,  If  v - 1.5,  the  in-sitc  SDF 
varies  as  f”^. 

To  determine  from  .9p(f ),  one  must  know  the  anisotropy  of  the 
irregularities  and  the  probe  velocity.  Moreover,  the  rms  electron 
density  Is  subject  to  the  same  ainbigulty  as  the  r.^8  phase.  However,  if 
we  assume  that  q^  * where  is  the  effective  length  of  data 

Interval  used  in  estimating  ?p(f),  then  Eq.  (7)  can  be  used  to  get  some 
ides  of  the  actual  density  perturbation  levels  involved. 

In  Figure  23  we  show  the  rms  electron  density  versus  for  data 

intervals  varying  logarithmically  from  l*km  to  1000  km.  Me  see  that 

even  for  the  200* km  layer  we  used  for  reference,  the  corresponding 

22 

glgahurts  level  of  10  gives  rise  to  rms  electron  density  utcursions 


L - aOQ  km 
H - »0  km 
X - 03  m (1  QHil 


of  loos  thiitt  10^^  ol/ni^  ovof  1 to  1000  ko  dots  lotorvolo.  ligaifleoatly 
•mallor  voluoo  reoult  If  tho  lrrotol**itioo  «ro  oxtondod  ovor  m lorgor 
roglon.  Thua,  oa  find  no  Inconalatoney  batwaan  our  aolntlllatlon 
obafiirvatlona  and  raportad  In*  at  to  Moauraawnta. 
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Appendix 

NOISE  COSSECTIONS  FOR  MEASURED  S^  SCINTILLATION  INDEX  VALUES 

If  we  apply  the  conventional  model  of  a alngal  plus  Independent 
additive  white  Gauseian  noise,  it  is  easily  shown  that  the  measured 
scintillation  index  has  the  form 

.2  1 1 

§4  “ (S4  + 2SNR‘^  + SNR-^)  /a  + SM”^>  (A-1 

wiiere  SNR  denotes  the  power  signal- to- noise  ratio.  It  is  interesting 
to  note  that  it  « 1,  then  §4  ” 1 irrespective  of  SNR, 

In  general,  however,  the  noise  contribution  causes  to  be  larger 
than  S4.  In  Figure  A-1  we  have  plotted  the  measured  scintillation  :|»ndiic 
against  S4  for  different  SNRs.  It  can  be  seen  that  when  < 0.1  ‘witli 
SNRa  less  than  30  dB,  a noise  correction  must  be  applied.  For  the 
Wideband  satellite  data,  the  SNR  is  typically  better  than  30  dB. 
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G. 
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